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Abstract 

Scattering amplitudes of gluons coupled with a pair of massive scalars, so-called massive 
scalar amplitudes, provide the simplest yet physically useful examples of massive amplitudes. 
In this paper we construct an S-matrix functional for the massive scalar amplitudes in a 
recently developed holonomy formalism in supertwistor space. From the S-matrix functional 
we derive ultra helicity violating (UHV), as well as next-to-UHV (NUHV), massive scalar 
amplitudes at tree level in a form that agrees with previously known results. We also obtain 
recursive expressions for non-UHV tree amplitudes in general. These results will open up a 
new avenue to the study of phenomenology in the spinor-helicity formalism. 



1 Introduction 



Recently there has been much progress in the computation of scattering amplitudes in four- 
dimensional massless gauge theories by use of the spinor-helicity formalism in twistor space. 
From technical and practical perspectives, most of the recent developments can be under- 
stood in a form of either the CSW rules [1] or the BCFW recursion relations [2, 3]. In order 
to apply these developments to phenomenological models, notably, in search of theories be- 
yond the standard model of particle physics, it is then natural to consider applications of the 
CSW/BCFW method to theories with massive particles. Indeed, such massive models were 
sought and investigated right after the proposals of these methods; for the case of the CSW 
rules, see [4, 5, 6] and for the BCFW relations, see [7]- [11]. For earlier works on electroweak 
phenomenology in terms of the spinor-helicity formalism, not exactly in a twistor framework, 
see, e.g., [12, 13, 14]. Some of more recent developments along these lines can also be found 
in [15]-[23]. 

Of these recent investigations the simplest massive models are presumably given by the 
scattering amplitudes of gluons coupled with massive scalars. These amplitudes, which we 
shall call massive scalar amplitudes from here on, are of direct relevance to one-loop calcu- 
lations in non-supersymmetric theories including QCD. Also, the massive scalar amplitudes 
are closely related to multigluon amplitudes with massive fermions, particularly quarks, by 
use of the supersymmetric Ward identities [24]. Thus a thorough and systematic understand- 
ing of the massive scalar amplitudes is crucial to build any phenomenological models in the 
spinor-helicity formalism. Some clues to such an understanding are already known in the 
literature. Particularly, Boels and Schwinn have obtained an analog of the CSW rules, the 
so-called massive CSW rules, for the massive scalar amplitudes [15, 16]. More recently, in 
[21] Kiermaier shows that the massive CSW rules correctly lead to the scattering amplitudes 
of a pair of massive scalars and an arbitrary number of positive-helicity gluons, the so-called 
ultra helicity violating (UHV) amplitudes, whose compact expressions have been derived pre- 
viously by BCFW-type recursion methods [7, 9, 11]. For the next-to-UHV (NUHV) massive 
scalar amplitudes, their CSW-type representations are essentially obtained by Elvang, Freed- 
man and Kiermaier (EFK) in the study of one-loop calculations for what is called one-minus 
amplitudes in QCD [22]. 

Motivated by these stimulating results, in the present paper, we consider construction of 
an S-matrix functional for the massive scalar amplitudes within the framework of a recently 
proposed holonomy formalism in twistor space [25]- [28]. There are a few good reasons to 
execute this study. First of all, in the holonomy formalism the CSW rules are implemented by 
a Wick-like contraction operator in a systematic functional language. This implementation 
is not limited to tree amplitudes; as demonstrated in [28], it can also be applied to one-loop 
amplitudes in J\f — 4 super Yang-Mills theory. Thus our primary concern is not the search of 
possible applications of the massive CSW rules to loop amplitudes. We would rather focus 
on the understanding of how the massive CSW rules are incorporated into the holonomy 
formalism at tree level, in expectation of how to obtain an insight into an utterly new mass 
generation mechanism. 

Secondly, a massive extension of the holonomy formalism is rather straightforward at least 
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from an algebraic perspective. As in the massless case, we need to define a massive holonomy 
operator so as to obtain an S-matrix functional for the massive scalar amplitudes. Practically, 
this can be carried out by making a massive extension of a bialgebraic comprehensive gauge 
field such that it satisfies the infinitesimal braid relations [29, 30]. As discussed in detail in 
section 3, it turns out that such an extension is indeed possible, which, in turn, algebraically 
guarantees the construction of the massive holonomy operator. 

Lastly, we notice that our construction is in accord with the recently studied on-shell con- 
structibility of massive amplitudes in general [18, 20]. In the holonomy formalism, physical 
information (i.e., helicity and a numbering index) is encoded in the creation operator of the 
involved particles. This principle should be held even for massive particles as there are no 
other ingredients for this role once a holonomy operator is defined. This implies that we can 
specify the polarization of a massive particle in a similar fashion to the case of helicity, i.e., 
we may also implement the polarization information into the massive creation operator by 
modifying Nair's prescription of superamplitudes [31]. Such a modification can naturally be 
made by an off-shell continuation of the null spinor momenta; notice that one can utilize the 
massive spinor-helicity formalism [14] to obtain an explicit form of massive spinor momenta. 
We shall confirm these interpretations in section 4 by presenting an S-matrix functional for 
the UHV massive scalar amplitudes. 

This paper is organized as follows. In section 2, we review the foundation of the holonomy 
formalism. Materials covered in this section are essential for later discussions. In section 
3, we show that the original massless holonomy operator can naturally be extended to a 
massive case from an algebraic point of view. We then define a massive holonomy operator 
for gluons and massive scalars. In section 4, we first consider off-shell continuation of Nair's 
superamplitude method and then briefly review the recent results of the massive CSW rules 
by Boels and Schwinn and their applications to the computation of the UHV massive scalar 
amplitudes by Kiermaier. We end this section by deriving an S-matrix functional for the 
UHV amplitudes in terms of the above obtained massive holonomy operator. In section 5, we 
extend the S-matrix functional to the NUHV amplitudes and confirm that our computation 
is in accord with the EFK result. We further discuss that the extended S-matrix also leads to 
recursive expressions for non-UHV massive scalar amplitudes in general. Lastly, we present 
concluding remarks. 

2 Review of holonomy formalism 

In this section, we review the foundation of the holonomy formalism introduced in [25] and 
developed in [26, 27, 28]. Materials covered here are indispensable for later discussions but 
the readers who are already familiar with the holonomy formalism may skip this reviewing 
section. 

Knizhnik-Zamolodchikov connections, twistor space and spinor momenta 

In the holonomy formalism, the holonomy operator refers to a holonomy of the so-called 
Knizhnik-Zamolodchikov (KZ) connection [29, 30]. The KZ connection in general is defined 
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by 

n = - Yl "a (2- 1 ) 

K i<*<i<™ 

where k is a non-zero constant, the so-called KZ parameter, and flij can be expressed as 

Qij = aS +) ® aJ- ) + 4~ } <g> a$ +) + 2aS 0) <g> af . (2.2) 
Here the operators and a l -°' ) (i = 1, 2, • • • , n) form the SL(2, C) algebra: 

[oJ +) ,a<->] = 2o{ 0) < 5y, [a<%«] = a< +) * y , [af\ a^} = -a f H <Sy (2.3) 

where Kronecker's deltas show that the non-zero commutators are obtained only when i — j. 
The remaining commutators, those expressed otherwise, all vanish. These operators act on 
a set of Fock spaces Vi which are characterized by the numbering indices %. In the holonomy 
formalism, the operators are identified with the creation operators of the i-th gluon 
with helicity ±. The physical Hilbert space of the holonomy formalism is then given by 
y®n = Vi <g) V2 <g) • • • <g) V^j. Vlij in (2.2) is a bialgebraic operator and its action on V® n can 
explicitly be written as 

^ 1 ® • • • ® 1 (8) Pi{I^) (g) 1 ® ■ ■ ■ (8) 1 (8) ® 1 ® • • • ® 1 (2.4) 

where J M (/i = 0, 1, 2) are elements of the SX(2, C) algebra, p denotes its representation and 
1 denotes the identity representation. 

The Ui/s in (2.1) are defined by the differential one- forms: 

Uij = d\og(zi - zj) = — - dZj (2.5) 

Zi Zj 

where the set of complex coordinates (i — 1,2, • • • ,n) are identified with local coordinates 
on CP 1 fibers of twistor space. In the holonomy formalism, these coordinates are related to 
the homogeneous coordinates of spinor momenta for the i-th gluon. The spinor momenta 
are parametrized in terms of null four-momenta for gluons. One of such parametrization is 
given by 

I) (2.6) 



>/Pi- Pi 



Pi ~Pi V Zi 



where A = 1, 2 and is a non-zero complex number, ctj e C — {0}. The null four-momentum 
Pi (p — 0, 1, 2, 3) satisfies the on-shell condition 

{Pi) 2 = ri^m = (P°) 2 - {p\f - (P 2 ) 2 - {Plf = (2.7) 
where we use the Minkowski signature (H ) for the metric i]^ u . 

Lorentz transformations of uf are given by uf — > guf where g G SL(2, C) denotes 
a (2 x 2)-matrix representation of SL(2,C). Scalar products of u A, s, which are invariant 
under the SL(2, C), are expressed as 

Ui ■ Uj = {uiUj) = e A Bufuf (2.8) 
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where €ab is the rank-2 Levi-Civita tensor. Similarly, we can define the scalar products of 
the complex-conjugate spinor momenta (A = 1,2) as 

tit ■ Uj = [UiUj] = ^ AB u iA UjB- (2.9) 

The null four-momenta are parametrized by the combination of the holomorphic spinor 
momenta uf and the antiholomorphic ones u i ^. In the spinor-helicity formalism, the four- 
dimensional Lorentz symmetry is therefore given by SL{2, C) x SL(2, C). 

In terms of the holomorphic spinor momenta, the logarithmic one-forms Uij in (2.5) can 
also be written as 

Uij = d\og(u iUj ) = -^y- (2.10) 

The physical configuration space of the holonomy formalism is given by = C n /S n 
where n is the number of gauge bosons, C™ represents a set of the Zi coordinates (i = 
1, 2, • • • , n) and S n denotes the rank-n symmetric group. The fundamental homotopy group 
is given by the braid group I1 1 (C^ A ' ) ) = B n . 

Infinitesimal braid relations and the integrability of KZ connection 

The integrability of the KZ connection, i.e., dQ — QAQ = 0, is guaranteed if fly satisfies 
the following conditions [29]: 

[flij^ki] = k, I are distinct), (2-11) 

[Vtij + Clj k , flik] = k are distinct). (2-12) 

These relations are known as the infinitesimal braid relations. The commutators of bialge- 
braic operators are generally defined by 

[oj ® bi, a,j <8> bj] = [oj, a,j] <S> h <8> bj + Oj <8> [h, dj] ® bj + a-i ® Oj <8> bj] (2-13) 

where a, and bi (i — 1, 2, • • • , n) denote a set of arbitrary operators. From (2.2) and (2.3), we 
find that the first relation (2.11) is obviously satisfied. One can also check that f^-'s satisfy 
the second relation (2.12). 

Comprehensive gauge one-forms for gluons 

Application of these mathematical results has lead to the holonomy formalism for gluon 
amplitudes. The physical operators of gluons are given by a^K Q^s are not appropriate 
to describe gluons since its action on the Hilbert space in (2.4) contains the action of af\ 
We need to modify O^-'s so that the operators af^ are treated somewhat unphysically, which 
leads us to introduce a "comprehensive" gauge one-form 

A = g Yl A i3 u l3 (2.14) 

l<i<j<n 

where g is a dimensionless coupling constant and A^ is defined as 

A i:j = <g> af + a [ r ] ® af . (2. 15) 



Notice that also satisfies the infinitesimal braid relations (2.11), (2.12); see [25] for 
details of its proof. As mentioned earlier, these relations guarantee the integrability of the 
"comprehensive" gauge field, i.e., 

DA = dA-AAA = -AAA = (2.16) 

where D denotes a covariant exterior derivative D = d — A. 

The coupling constant g is related to the KZ parameter k by g = K For an SU (N) gauge 
theory, this can be given by 

g = - = — *— . (2.17) 
y k 1 + N v ; 

Definition of the holonomy operator for A 

The integrability of the comprehensive gauge one-form A allows us to define a holonomy 
of A. The holonomy operator of A is defined by 



® ( rI( u ) = Tr^Pexp 



<f A A A A ■ ■ ■ A A 

m>2 J 7 " m ' 



(2.18) 



where 7 represents a closed path on = C n /S n along which the integral is evaluated 
and R denotes the representation of the gauge group. The color degree of freedom can be 
attached to the physical operators a- ± ' ) in (2.15) as 

aS ±) =t-aS ±) - (2.19) 

where t Ci 's are the generators of the SU(N) gauge group in the /^-representation. The 
symbol P denotes an ordering of the numbering indices. The meaning of the action of P on 
the exponent of (2.18) can explicitly be written as 



m>2 J ~1 

E f A 12 A 23 ■ ■ ■ A ml u 12 A co 23 

m>2 J ~< 

E ^1 E (-l) fcl+fca+ - +fc » a? ® 4 h2) ® • • • ® j> u 12 A • • • A u . 



• • • A A 

A • • • A u ml 



m>2 (hi,h2,—,hm) 



(2.20) 



where hi = ± = ±1 (i — 1, 2, • • • , m) denotes the helicity of the i-th gluon. In deriving the 
above expression, we use the relations 



[^12,^23] 

[[Ai 2 ,A 23 ], A u ] 



4 +) ® 


4 +) ® 


a {0) - 
a 3 


4 +) ® 


4 ) ® 


a (0) 

°3 








® 4 +) 


®4 0) 


-4" } 


® 4~ } 


®4 0) 






4 +) ® 


4 +) ® 


4 +) ® 


4 0) - 


4 +) ® 


4 +) <g> 


4 _) ® 




-4 +) 


<S> 4 _) 


® 4 +) 


®4 0) 


+ 4 +) 


® 4~ } 


® 4 _) 


®4 0) 


+ 4" } 


® 4 +) 


<g> 4 +) 


®4 0) 


-4- } 


® 4 +) 


<g) a 3 _) 


®4 0) 


+ 4") 


<g> 4 _) 


® 4 +) 


®4 0) 


-4- } 


<g> 4 _) 


® 4~ } 


®4 0) 



(2.21) 



(2.22) 
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and their generalization. In the expression (2.20), we also define (gm^ 



4 hm) ®4 0) 



as 



<g> a, 



(km) 



4 0) - ^[ai 0) ,aW]®^ 
= ±iaS ±} (g) af 2) (g) • • 



<g> a. 



(/in 



(g a 



(2.23) 



where we implicitly use an antisymmetric property for the numbering indices (1, 2, • • • , m). 

The trace Tr# i7 in the definition (2.18) represents a combination of the usual color trace 
Tir over t Ci 's and the so-called braid trace Tr 7 over braid generators. The braid trace is 



realized by a sum over permutations of the numbering indices; see [26] for details of this 
point. Thus the braid trace Tr 7 over the exponent of (2.18) can be expressed as 

Ar m i u la2 A u a2U3 A • • • A u Uml (2.24) 



m>2<reS m _i T 



where the summation of <S m _i is taken over the permutations of the elements {2, 3, • • • , m}, 

/ 2 3 • • • m \ 

with the permutations labeled by a = 

v (T 2 (73 • • • a m J 

The holonomy operator in supertwistor space 

In the holonomy formalism, an S-matrix functional for gluon amplitudes is described by 
a holonomy operator in supertwistor space. The supersymmetrized holonomy operator is 
defined by 



ei#(u;M) = Tr^Pexp 



V / A A A A ■ ■ ■ A A 



where the bialgebraic operator A^ in (2.15) is now expressed as 



hi 



(2tt) 3 2^ 



//, • r/ // . (f 2 II ; 



Ui ■ dui d Ui 
2ni (2tt) 2 2rri (2vr) 2 



(2.25) 

(2.26) 
(2.27) 
(2.28) 



dfi(pi) is called the Nair measure for the null momentum p^. a[ hl \x, 6>)'s are physical operators 
that are defined in a four-dimensional Af = 4 chiral superspace (x, 6) where x AA denote 
coordinates of four-dimensional spacetime and 8 A (A = 1,2; a = 1,2,3,4) denote their 
chiral superpartners with M = 4 extended supersymmetry. These coordinates emerges from 
homogeneous coordinates of the supertwistor space CP 3 ' 4 , represented by (u A , v A , £ a ), that 
satisfy the so-called supertwistor conditions 



V A 



X AA U 



e a A u A 



(2.29) 
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The physical operators af 1 ^ are relevant to creations of gluons and their superpartners, 
having the helicity hi = (0, ±|, ±1). Explicitly, these supermultiplets can be expressed as 



4 +) te) 


= 4 +) , i 








. 1 pccpV JO) 











(2.30) 



which are consistent with the definition of the helicity operator 



(2.31) 



Use of the supermultiplets (2.30) enables us to define gluon operators without introducing 
the conventional polarization/helicity vectors. This method is known as Nair's prescription 
of superamplitudes [31]. 

An S-matrix functional for gluon amplitudes 

In terms of the supersymmetric holonomy operator (2.25), an S-matrix functional for 
gluon amplitudes can be constructed as 



where 



MHV 



(h)c 



W {A \x) 



= exp 



cxp 



exp 



(h)c 



MHV 



(h)c 



- 2 Jd 4 xd 8 e eS#( u; x,e) 



-J dfi(q) 



(2.32) 



(2.33) 



5a ( P +) 



5a 



(-) 



-iq-(x-y) 



d 4 q 



(2tt)V UoJ- 



(+) 



5a 



(-) 



y-tx 
-iq-(x-y) 



(2.34) 



y-tx 



Note that we take the limit y — > x, keeping the time ordering x° > y° or x° — y° — > 0+, 
at the end of calculation. The CSW rules are realized, in a functional language, by the 
incorporation of the Wick-like contraction operator W^ A \x) into the S-matrix functional 

jr(A) Uh)c] 



. In (2.34), q denotes a momentum transfer which is generally off-shell and p 

(2.35) 



denotes its on-shell partner. The two are related by 

q» = Pv + wr]^ 



where 77^ is a reference null-vector, satisfying rj 2 = and w is a real number. Since both 77^ 
and w can arbitrarily be chosen, we can fix the scaling freedom for either 77^ or w. 
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In terms of the S-matrix functional T a^ c , general n-point N fe MHV gluon amplitudes 
(k — 0, 1, 2, • • • n — 4) are generated as 



6 6 .jr(A) 



Sa{ hl)ci 5a ( 2 h2)c2 5a ( n hn)Cn 



4Zuw nhn \x) (2-36) 



a C0 c =0 



where a^ c denotes a generic expression for the gluon creation operators af~^ c% (with hi = ±, 
i = 1, 2, • • • , n), which are treated as source functions in the above. Notice that the expression 
(2.36) is not limited to the case of tree amplitudes. As shown in [28], the expression is also 
applicable to one-loop amplitudes and, from a functional perspective, it would and should be 
valid through higher loop levels. 

In practical calculations, we need to use two key relations. One is the normalization of 
the spinor momenta 

/ d{u lU2 ) A d(u 2 u 3 ) A • • ■ A d(w m «i) = 2 m+1 (2.37) 
and the other is the non-vanishing Grassmann integral over 6*'s: 



;a =9 a A 

n A % 



u r u s y . (2.38) 



The latter relation guarantees that gluon amplitudes vanish unless the helicity configuration 
can be factorized into the MHV helicity configurations. Together with use of the contraction 
operator (2.34), the CSW rules are thus automatically satisfied by the Grassmann integral 
(2.38). 

Lastly, to clarify the notations above, we present the tree-level MHV amplitudes, the 
simplest form of the gluon amplitudes, in the x-space representation [25]: 

^MHV(0) + \x) = »4mhv(o) ( x ) = IT / dfJ>(Pi) •^MHV(0)(' U ' ^) ' (2.39) 

i=l 

<mvi)(^) = ^-W^ (4) (i>) A^ 0) (u), (2.40) 
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Of course, there exists a lot of complexity in the generalization of these forms to non- 
MHV and higher-loop amplitudes but from the above functional expression (2.36) we can in 
principle write down the x-space N fc MHV gluon amplitudes as [28]: 

A (l hl 2 h2 -n hn ) , , _ (l hl 2 h2 -n hn ) , , (l hl 2 h2 -n h „), , .(l hl 2 h2 -n hn ) , , , - 

N k MHV \ X J ~ N k MHV(0) \ X > " r ,/H N k MHV(l) v X > " r ,/i N k MHV(2) \ x > ~T~ " ' V^-^J 

where A^^^y^ n \x) denotes the n-point L-loop N fc MHV gluon amplitude and hi — ± 
denotes the helicity of the i-th gluon, with the total number of negative helicities being k + 2 
(fc = 0,l,2,...,n-4). 
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These are basic results of gluon amplitudes in the holonomy formalism. Since we consider 
a purely gluonic theory, the helicity index is specified by hi = (+,—), rather than the 
supersymmetric version hi = (0, ±±, ±), as shown in (2.36) and (2.42). If we include massive 
scalars, however, we need to incorporate hi = ingredients due to the definition of the 
helicity operator (2.31). Consequently, it is inevitable to modify the purely gluonic S-matrix 
functional (2.32). In order to implement such a modification, in the next section we first 
consider a massive extension of the holonomy operator. 



3 The holonomy operator of gluons and massive scalars 

In this section, we consider incorporation of massive operators into the holonomy operator 
from an algebraic perspective. The aim of this section is to construct a holonomy operator 
that is relevant to the massive scalar amplitudes, i.e., the amplitudes of gluons coupled with 
massive scalar particles. 

Massive extension of comprehensive gauge fields 

To begin with, we consider a massive extension of the comprehensive gauge field A in 
(2.15). We find that the most natural extension can be made by 1 

B= £ SyWy (3.1) 
l<i<j<n 

where the "massive" bialgebraic operator By is given by 

B^ = g (oS +) ® af + oS-> ® af) + af ® of . (3.2) 



As before, Uij is the logarithmic one-form in (2.10) and g denotes the dimensionless gauge 
coupling constant. 

From the definition (3.2), one can easily check that B^ satisfy the infinitesimal braid 
relations: 

[Bij,B k i] = k, I are distinct), (3.3) 

[By + Bj k , B ik ] = k are distinct). (3.4) 

The first relation (3.3) is trivial from (2.3) and (2.13). The second part can also be checked 
by 

[Bij,B ik ] = <? 2 [<4 +) ®af ® af]+ g[af ] ® af ,af ® af] 
+ 9 2 k H ® af, a« 4 0) ] + g[a^ } ® af,af ® 4 0) ] 

(°) „(0) „(+) ^ „(0)i , „rA0) ^ (0) (-) ^ (Oh 



+g[ar®af,ar > ®ar]+9[ar®af,a 



k 



= (3.5) 



1 The choice of ft in (2.1) also seems reasonable at first glance since, as discussed in the previous section, 
it satisfies the infinitesimal braid relations. But calculations of [Sl 12 , ^23], [[^12, ^23], ^34], etc., indicate that 
the holonomy operator of fi leads to unwanted prefactors. 
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and the trivial relation [Bjk,Bik] = 0, with the indices k being distinct. 

Definition of a holonomy operator for B: a first look 

Since the infinitesimal braid relations are satisfied, we can naively define a holonomy 
operator of B as 



e5gH=TY R , 7 Pexp 



j B A B A ■ ■ ■ A B 

r>2 J T ^ ? ' 



(3.6) 



As in (2.22) and (2.23), an explicit expansion of the physical operators a\ h ^ (hi = ±, 0) 
in the integrand can be deduced from the commutation relations 

[B 12 ,B 23 ] = g 2 (4 +) ® a { 2 +) ® af - a[ +) <g> a ( f ] ® af 

+aS- ) ®4 +) ®4 0) -aS- ) ®4- ) ®4 0) ) 

+g (aS 0) ® 4 +) ® af - ® 4 _) ® af) , (3.7) 
[[5i 2 , 5 23 ] , 5 3 4] = # 3 (4 +) ® 4 +) ® 4 +) ® af - 4 +) ® 4 +) ® ® af 

- 4 +) ® ® 4 +) ® ai 0) + 4 +) ® 4 _) ® 4~ } ® 4° } 
+ 4~ } ® 4 +) ® 4 +) ® 4° } - 4 _) ® 4 +) ® 4 _) ® 4° } 

- 4~ ) ® 4 _) ® 4 +) ® 4 0) + 4 _) ® 4~ } ® 4 _) ® 4 0) ) 



V (4 



" ,] ® 4 +) ® 4 +) ® af - ai 0) ® 4 +) ® a<-> ® af 

f } <g) 4 _) ® 4 +) ® 4 0) + 4° } ® 4 _) ® 4 _) ® 4° } ) (3.8) 



— a 

and their generalization. Terms in the leading order of g are the same as the massless case. 
Thus, by use of the definition (2.23), these terms lead to the original massless holonomy 
operator, reducing the holonomy of B to that of A. 

The rest of the terms, those with a^'s, would correspond to correlators of the interaction 
among gluons and a pair massive scalars. Note that, as discussed in (2.30) and (2.31), a 
creation operator of a scalar or spin-0 particle is described by af^ in the holonomy formalism. 
If we apply the definition (2.23), however, these would-be massive terms vanish and we can 
not construct a massive holonomy operator out of (3.6). This problem can be remedied by: 

1. considering an open path integral so that a pair of the operators (4°\ a ffl) survives in 
the integrand of (3.6); and 

2. splitting the numbering indices into those of gluons and massive scalars when we take 
a braid trace or a sum over the permutations of indices. 

For this purpose, we first fix the indices of massive scalars to 1 and n, being in accord with 
the expressions (3.7), (3.8). We then identify 4° } and 4°) as physical operators for a pair 
of complex massive scalar particles 0i and </> n , respectively 2 . We consider that gluons and 

2 We here follow the convention to use complex particles. As we shall see later, no significant differences 
arise between real and complex massive scalars in our formalism. Use of complex scalars is simply more 
suitable for the extension to amplitudes of gluons and fermions. 
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massive scalars are both in the ^-representation of the gauge group, SU(N) xU(l) = U(N), 
otherwise we can not properly define couplings between them. Notice that, as in the case 
of scalar propagators that appears in the CSW rules, we can assign U(l) color degrees 
of freedom to the scalar particles so that the single trace structure of the full amplitudes 
preserves. 

Consequently, the braid trace in (3.6) should be taken over the numbering elements 
{(7 2 , <7 3 , • • • , oy_i, r r } = {2, 3, • • • , r}, satisfying the P ordering 

a 2 < a 3 < ■ ■ ■ < ay_i . (3.9) 

The braid trace can then be represented by a "homogenous" sum 

E = E E (3.io) 

{<t,t} r r =2 <reS r _2 

where r = 3, 4, • • • , n. 

As discussed in [27] (see section 3), a product of iterated integrals over the logarithmic 
one-forms w^-'s can be expanded, using the homogeneous sum, as 

^ iu lff2 Aw w A---Aw ffr _ 1Tr Au Tr i = / w 12 A • • • A w r _ ir / u rl (3.11) 

where 7i r and 7 r i denote open paths on a physical configuration space of interest, satisfying 
7 — 7ir7ri- Notice that we split the numbering indices into <7j (i — 2, 3, • • • , r — 1) and r r , 
respectively corresponding to the elements of gluons and a pair of massive scalars. In this 
labeling, the closed path 7 can be denoted as 7 = 7u | t * The physical configuration is now 
given by that of (n — 2) gluons and 2 distinct massive scalars, i.e., 

C (B) = ^_ C 2 = C n /Sn _ 2 (3 . 12) 

as opposed to the pure gluonic case 

C (A) = Q n /S n . In the present massive case, any physical 
observables should be symmetric under transpositions of (n — 2) gluons. This is consistent 
with the appearance of the sum over a G 5 r _ 2 in (3.10). The quantum Hilbert space, on the 
other hand, remains the same as the massless case, V® n — V\ ® V 2 <8> • • • <S> V n as discussed 
below (2.3). 

Definition of a holonomy operator for B: a refined version 

Now that we have specified the physical configuration on which the massive holonomy 
operator Q^(u) is defined and the quantum Hilbert space on which 6^(w) acts, we are at 

the stage of deriving a well-defined version of 6^(m) to replace the naive guess form (3.6). 

From the above arguments, we find that an analog of the expansion (2.24) can be ex- 
pressed as 

00 . 

Tr 7 P f B A ■ ■ ■ A B 

12 



— E E f B la2 B a2(T3 ■ ■ ■ B ar _ lTr B Trl co la2 A uJa 2 a 3 A • • • A LU ar _ lTr A to Tr i 
= E / B 12 B 2 3 • • • B r _ lr u>i2 A • • • A ov_i r / _B rrl w rl 

r>3 



(3.13) 



where we treat B^s as coefficients of the logarithmic one-forms. Since the pure gluonic part 
is excluded from the physical configuration space (3.12), the above integral leads to operators 
involving gluons coupled with a pair of massive scalars (a^\ a^). Namely, we have 



Tr 7 P f BA- - -AB 

E E f' 2 

r>3 (h 2 ,h 3 ,---,h r -i) 



h2hs-h r . la (0)^ a (h 2 ) 



x 



/ 



U\2 A'-'Aw. 



■-lr / 

•>7rl 



u; r i 



(0) „ (h 2 ) 

= E E g r - 2 (^-l^ h 3-hr-l^l 



r>3 (h2,hz,-,h r -i) 



(12)(23)---(r- lr)(rl) 



(3.14) 



where /ij = ± = ±1 (i = 2, 3, • • • , r — 1) and we use the normalization of the spinor momenta 
along an open path 



/ d(uiu 2 ) A d(u2Us) A • • • A d(u r -iu r ) = 2 r 1 , 

/ d(u r Ux) = 1 . 

Jjrl 



(3.15) 



These are open-path analogs of the closed-path normalization given in (2.37). 

To summarize, we can define the holonomy operator of gluons and massive scalars as 



e S( u ) = ex P 



T £ ^-2(_ 1 )feft3-^-l Xr ^C 2 ^3... t C r _ 1 ^ 

r>3 (h2,h 3 ,-,h r -i) 



X 



aS 0) ® a^ 2)c2 



a 



(h r -i)c r 



®a£ > 



(12)(23)---(r- lr)(rl) 



(3.16) 



where we make the color factor explicit. Notice that the braid trace, or a sum over permu- 
tations of gluons, is not apparent in this form but it is already taken account of in splitting 
the original closed path 7 into two open paths 7i r and 7^1 . Thus the braid trace is implicitly 
realized by the homogeneous sum (3.10) with use of the relation (3.11). This form is different 
from the conventional color decomposition of the massive scalar amplitudes where the sum 
over permutations is explicit as in the pure gluonic case; see, for example, [24] and references 
therein. 
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4 An S-matrix functional for UHV tree amplitudes 



As mentioned earlier, in the holonomy formalism all the physical information should be 
encoded in the creation operators, i.e., in af^'s for gluons and in (a^\a^) for massive 
scalars. In the case of gluons, the helicity information is implemented by supersymmetriza- 
tion of the underlying twistor space. As discussed in section 2, this is implemented by Nair's 
superamplitude method [31]. In this section, we first consider off-shell continuation of this 
method. We then briefly review the recent results of the massive CSW rules [15, 16] and 
their applications to the computation of the so-called ultra helicity violating (UHV) ampli- 
tudes, i.e., the scattering amplitudes of a pair of massive scalars and an arbitrary number 
of positive- helicity gluons, at tree level [21]. To the end of this section, we shall present 
an S-matrix functional for the UHV tree amplitudes by introducing a Wick-like contraction 
operator involving the massive operators. 

Off-shell continuation of Nair's superamplitude method 

To begin with, we rewrite the off-shell parametrization of a four-momentum (2.35) as 

where denotes a massive four-momentum with mass m and p M denotes its on-shell partner. 
rf is a reference null- vector. In terms of spinor momenta, the null momentum is expressed as 
p AA = (<Jfj) AA p^, with A and A taking values of (1,2). here is given by = (1, a A where 
<7j (i = 1,2,3) and 1 are the Pauli matrices and the (2 x 2) identity matrix, respectively. 
Using the parametrization (4.1), we can then define off-shell continuation of the null spinor 
momenta as [24] 

u A — )■ u A = u A + -^-r] A , (4.2) 

(urj) 

u A — > U A = u A + ^-n A (4-3) 

[ufj] 

where r/ A is a reference null spinor and f\ A is its complex conjugate. 

Since the reference null-vector i] AA = i] A f] A can be chosen arbitrarily, it is defined on 
a distinct twistor space, decoupled from the original one that has been parametrized by 
the spinors (u A ,Vj) satisfying the condition v A = x^ A u A . This interpretation of i] AA is in 
accord with the definitions (4.1)-(4.3). In order to construct a massive model in the spinor- 
helicity formalism, however, naive substitution of w A 's by u A, s does not work out well. For 
example, one can consider that an off-shell continuation of the the projected Grassmann 
variable £ Q = Q\u A in (2.29) is given by \ a = 9^u A . Use of £ a in the expressions of Nair's 
superamplitude method (2.30) leads to vanishing UHV amplitudes due to the Grassmann 
integral (2.38). But this is contradictory because, as reviewed below, the UHV amplitudes 
are non- vanishing in general. 

Simple use of (4.2) and (4.3) therefore does not lead to massive extensions in the spinor- 
formalism. In fact, one should rather think of two distinct sets of twistor variables (u A ,v^) 
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and 



{W A , TT A ) 



where w A = -r^r\ A and 1: A = x AA W 



spinor basis spanned by [20] 



(U7j) 



In other words, we should use a two- 



u 



m 



(urj) 



m 



u 



(4.4) 



to describe holomorphic and antiholomorphic massive quantities, respectively. 



Notice that the four-dimensional spacetime x AA emerges from each of the twistor vari- 
ables. This feature should be preserved after supersymmetrization of the underlying twistor 
spaces. Namely, in addition to the original supertwistor variables {u A , v A , £ a ), we need to 
introduce new supertwistor variables {w A , n A , ( a ) such that the supertwistor conditions 



7I\ 



x AA w 



X 



AA 



m 

(urj) 



rj 



Q\w J 



0<: 



m 



(urj) 



rj 



(4.5) 



are satisfied (a = 1,2,3,4). The emergent chiral superspace, i.e., the four- dimensional 
spacetime x AA and its chiral superpartner 8 A , is identical for either the original or the new 
supertwistor spaces. This is explicitly presented in (2.29) and (4.5). 

We now consider off-shell continuation of Nair's superamplitude method. Based on the 
above arguments, this can be implemented by modifying the operators (2.30) for massive 
scalars in terms of £ Q 's and £ a 's. For this purpose, we take account of the conditions 
that (a) the UHV tree amplitudes are non-vanishing and (b) the massive scalar operators 
have 2 degrees of homogeneity in -u's. The latter condition is in accord with the helicity 
operator (2.31). Regarding the former, we shall show an explicit form of the UHV tree 
amplitudes later; see (4.28). Using the Grassmann integral (2.38), we then find that the 



massive scalar operators can uniquely be determined as c4^g(£i,C 



M 

M 



= \ifi-afl p where 

Y2 e a/3js£,i (f af ' 1 ■ In other words, we can define an off-shell continuation of the operator 



a) Uj (&) in (2.30) as 



a 



(0) 



(6,0) = a 



(o) 



(4.6) 



where we shall specify the numbering index to % — l,n for massive scalars. 
On the other hand, the gluon operators remain the same as in (2.30), i.e., 



4 +) 

.(-) 



,(+) 



(6) = 
(6) = £&f£<4 _) 



(4.7) 
(4.8) 



where % = 2, 3, • • • , n — 1 (with n = 3, 4, • • •). The gluonic part of the massive holonomy 
operator (3.16) is then automatically obtained by use of the Af = 4 chiral superspace rep- 
resentation (2.27) with the on-shell Nair measure (2.28). For the massive scalars, the same 
superspace representation can be obtained by use of off-shell continuation of the Nair mea- 
sure d/j,(pi). (Although we shall not use the off-shell Nair measure explicitly in the present 
paper, interested reader may refer to details of the off-shell Nair measure in [28].) The chiral 
superspace representation of the massive operators can then be expressed as 



,(°) 



(x,9) = j 'dfifa) aS 0) (6,Ci) 



=0%wf 



(4.9) 
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m 



Pt = rf + 55 ^, (4-10) 



[UiTJi) 

We can use the expressions (2.27) and (4-9) for gluons and massive scalars, respectively, 
to construct a super symmetric version of the massive holonomy operator. Namely, we can 

( B} ( B) 

obtain the super symmetric massive holonomy operator Q R ^(u;x, 9) out ofQ R ^(u) in (3.16) 

with replacements of {a^, a!f^} by {af^^x, 9), af\x, 9)} where i — 2, 3, • • • , r — 1 and j = 
l,r. 

Review of the massive CSW rules 

In the following, we briefly review the massive CSW rules of Boels and Schwinn [15, 16]. 
These are an analog of the original CSW rules for gluons gf and massive complex scalars 4>u 
(pi. As in the original case, the massive CSW rules give prescription for amplitudes in terms 
of vertices connected by massless and massive scalar propagators, 

D 9+ Af) = J^ D^f) = w ^— 2 (4.12) 

for positive and negative-helicity gluons and a pair of massive scalars, respectively. Up to 
constant factors, the involving vertices are expressed as 

VMmigtgt ■ ■ ■ 9i+i9z gf+i ■ • • 9t-i9j at+i • • • &T) = (i 2 )(23)---(ni) ' ^' 13 ^ 

(nl) 2 

W^ 2 + -^A) = ™ 2 (12)(23) ... (wl) (4-14) 

where Vmhv is a purely gluonic MHV vertices, with its form being the same as the orig- 
inal CSW rules. Peculiarity of the massive CSW rules lies in the form of Vuhv which is 
proportional to m 2 . There also exist non-UHV type vertices involving the massive CSW 
rules 



^NUH V (01#2~ • ' • 9i+l9i ft+i ■ ■ ■ 9n-M 
(li) 2 (in) 2 



(12)(23).-.(nl) 



(4.15) 



^UHV 2 (^92 ■ ■ ■ at+lMi+igt+2 ■ ■ ■ #n-l ( 



1 / (li) 2 (t + ln) 2 (li)(t + ln)(l» + l)(m) \ 

2 V(12)(23)---(nl) (12)(23)---(nl) J K ' ' 



which are not of direct relevance to the calculations of the UHV massive scalar amplitudes. 
In the above expressions the spinor momenta corresponding to the massive scalars, i.e., 
and are given by the on-shell partners of the actual massive spinor momenta u± and 
respectively. Notice that these spinor momenta are related to each other by (4.2). 



The massive CSW rules have been proposed by use of two Lagrangian-based methods. 
One is to use a canonical transformation in the light-cone gauge [32, 33, 34] and the other 
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is to use an action constructed in twistor space [35, 36]. In either approach, one starts from 
the ordinary Lagrangian for gluons and massive scalars 

1 



C = -— F ^ FCiiV + (D^) C (D^) C - ™W (4-17) 



c 

\1V 



where _D M = d^ + A^ (A^ = —it c A c ^) is the covariant derivative and F M „ = [D^, D v \ = —it c F^ 
is the field strength for gluons. Following the notation (2.19), we here denote the color 
factor by t c 's. One then carries out field redefinitions of and <fi such that non-MHV type 
couplings between a pair of massive scalar fields and gluons are eliminated. Together with 
supersymmetry arguments, this enables one to obtain the above forms of vertices [16]. 

Since the massive CSW rules are based on the Lagrangian formalism, they are not nec- 
essarily compatible with our holonomy formalism where we do/can not introduce massive 
potentials. In fact, using our parametrization (4.6) and the Grassmann integral (2.38), we 
can readily find that the vertices (4.15) and (4.16) vanish, i.e., 

VNUHv(01^-"^l^ r ^l"-^n-10n) = °, ( 4 - i8 ) 

V V uv*{(l>igt ■■■9t-i ( Pi < t ) i+i9t+2---9n-i ( Pn) = (4.19) 

where we omit the suffix (BS) to distinguish the vertices from those of the massive CSW 
rules. These results seem to contradict each other. In fact, however, although it is not 
well-recognized in the literature, there are no explicit derivations of the non-UHV type 
vertices V^uhv and V^y 2 , as clearly stated in [16] (see at the end of subsection 3.3). The 
NUHV vertex V^uhv) together with the UHV vertex Vuhv> does lead to four- and five-point 
NUHV massive scalar amplitudes [15, 16] but there are also possibilities that different NUHV 
vertices would lead to correct NUHV tree amplitudes because these amplitudes are generally 
dependent upon reference spinors which we can arbitrarily choose. We shall come back this 
point and consider such possibilities in the next section; see discussions below (5.8). 

As far as the UHV amplitudes are concerned, this apparent discrepancy goes away. To 
see this assertion, we now present a functional derivation of the UHV vertex (4.14) in terms 
of the massive holonomy operator (3.16). 

Choice of reference spinors and functional derivation of the UHV vertex 

We first fix the reference null- vector corresponding to the pair of massive scalar particles 

by 

rf=Pn, <=Pi- (4-20) 



This means that we have 



m 2 



Pi = Pi + tt ( zrt=l% + wrfl, (4-21) 

2(Pl • Pn) 



m 



PZ = Pi + 2 (p n . Pl ) P " =P " + WP " ' (422) 

satisfying p\ = p 2 n = rn 2 and m = 2 (pf- P y Since both pi and p% are massive, we can 
parametrized them as (4.21) and (4.22) in a suitable reference frame. This parametrization 
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is qualitatively different from off-shell prescription for virtual gluons where we set all reference 
null-vectors identical. 



Fixing the reference spinors as such, we now derive the UHV vertex (4.14) from the 
supersymmetric version of the massive holonomy operator (3.16) by a functional method. 
As in the MHV amplitudes, we introduce a generating functional 



-(vertex) 



= exp 



•> UHV 

Then the UHV vertex can be generated as 

5 5 5 5 



(4.23) 



5a\ 



(o) 



5ai +) 



5a 



(+) 



(vertex) 



5a { +\ 



5a 



(0) ^uhv 



a(±) c =a(°)=0 



(7 + + i \ r n—l 

J i=2 



(4.24) 



71-1 



-19" 



,n-2 



(27r) 4 ^)(p 1 + E^+Pn) V^u), 



i=2 



U, 



(<t>l<l>n) 



UHV 



(«) 



m 2 (nl) 2 



(12)(23)---(n- lra)(nl) 
= Tr(t c H C3 ■ ■ ■ t c ^) U UHV (<M 2 + • • • gtM 



where we use the Grassmann integral 



(In) 3 / 



m 



= (ln) J 



2 



m 



]2 Q 4 



;Vn0B 



(uiu n )(u n ui) J 



J d 2 6uie A ufe B = m 2 (nlf 



■ (u n ui) 



(4.25) 



(4.26) 



(4.27) 



This Grassmann integral guarantees that only the UHV-type vertices survive upon the eval- 
uation of functional derivatives in (4.24). This also automatically leads to the vanishing of 
non-UHV vertices (4.18) and (4.19), rather than (4.15) and (4.16). 

Another interesting feature in the expressions (4.24)-(4.27) is that there arise no sums 
over permutations of the numbering indices, contrary to the case of MHV gluon amplitudes 
(2.39)-(2.41). This implies that the number of terms to describe the UHV massive scalar 
amplitudes drastically decreases from that of the MHV gluon amplitudes. However, such 
a reduction does not occur in the massive scalar amplitudes. This is due to the fact that 
we can construct the UHV amplitudes by connecting the UHV vertices with as-many-as- 
possible massive propagators D^p) in (4.12). Notice that the number of propagators or 
vertices is independent of the gluon helicity configurations in the present case, while in the 
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gluon amplitudes the number of massless propagators D g + g -(p 2 ) in (4.12) is fixed by the 
helicity configurations or by the number of negative-helicity gluons. 

We can then express the UHV amplitudes by a UHV vertex expansion. As we shall review 
in a moment, indeed, such an expansion is explicitly realized in Kiermaier's expression for 
the UHV amplitudes [21]. Once the UHV amplitudes are constructed in this way, extension 
to next-to-UHV (NUHV) amplitudes which contains one negative-helicity gluon in addition 
to the UHV configuration, is straightforward by application of the original CSW rules or 
the MHV rules to the gluonic part of the amplitudes. Generalization to N fc UHV amplitudes 
{k — 1, 2, • • • , n — 3) can be carried out in the same manner. We can therefore construct the 
scattering amplitudes of an arbitrary number of gluons in any helicity configurations and a 
pair of complex massive scalars. We may call this construction the "UHV rules" for massive 
scalar amplitudes in analogy to the MHV rules for gluon amplitudes. We shall consider the 
non-UHV type constructions in the next section. 

The UHV expansion: Kiermaier's result for the UHV tree amplitudes 

Recently Kiermaier shows that the UHV tree amplitudes can be obtained by use of 
the massive CSW rules, more precisely, by use of the UHV vertex (4.14) and the massive 
propagator (4.12). The resultant expression is given by [21] 



-m 



^uhv(o) - (i2)(23)-..(n-ln) 



n-2 

n 

j=2 



1 - 



im 2 \J)(jj + l)(J\ 



{P 2 -m 2 )(jJ)(Jj + l)\ 



n (4.28) 



where (l|n) = (In) = eabu^u^ and 

P/ = Pi + P2 + P3 H 1- Pj ■ 

As before, the on-shell partner of Pj is defined as 



(4.29) 



2(pj • Vj) 



(4.30) 



where r/j denotes a reference null-vector. The corresponding spinor momenta (uj,Uj) are 
then defined by 

pj A = ujuj. (4.31) 



While the form (4.28) is probably the most concise expression of the UHV tree amplitude, 
for the clarification of the above mentioned "UHV rules," we now rewrite it as follows: 



i (0192 -9„-l4>n) 

l UHV(0) 



m 2 {n\) 



(12)(23)-..(nl) 



(nl) 



(4.32) 



where 



(nl) 



{nl) + j2 iJ1)im2{33 + 1) {nJ) 



i=2 



(jJ) P]-m 2 (Jj + 1) 
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+ E 



(II) im 2 (ii + 1) (JI) im 2 (jj + l) (nJ) 



+ E 

2<i<j<k<n-l 



+ 



(II) im 2 (ii + 1) (J/) im 2 (jj + l) 
(U) P 2 -m 2 (Ii + l)(jJ) P 2 -m 2 



x 



(KJ) im 2 (kk + l) (nK) 
(Jj + l)(kK) Pl-m 2 (Kk + 1) 



(4.33) 



Note that the factor im 2 in the numerators should be regarded as im 2 = \J^-\m 2 . The 
uppercase letters I, J,K, - ■ ■ play the same role as the J in (4.29)-(4.31). More explicitly we 
can expand the UHV tree amplitudes as 



A 



UHV(O) 

m 2 (nl) 2 
(12)(23)---(nl) 

m 2 (Jl) 2 i m 2 (nJ) 2 

+ ^ (12) (23) ■■■{jj- 1)(JV)(J1) p] - m 2 (J J + l)(j + 1 J + 2) • • • (n - 1 n)(nJ) 



+ E 

2<«<i<n-l 



m 2 (/l) ; 



m 2 (J/) 5 



+ E 

2<i<j<k<n-l 



12) • • • pa _ m 2 ( 7i + + u + 2) • • • (j - lj)(jJ)(JI) 

i m 2 (nJ) 2 
X Pj - m 2 (Jj + l)(j + 1 j + 2) • • • (n -ln)(nJ) 
m 2 (Il) 2 % m 2 (JI) 2 



(12) • • • pa _ m 2 ( 7i + i)(i + u + 2 ) . . . (j - 1 3 )( 3 J)(JI) 

% m 2 (KJ) 2 



x 



x 



Pj _ m2 (j j + + i J + 2 ) ■ ■ ■ (k - 1 k)(kK)(KJ) 



m 2 (nKf 



z(Kk + l)(k + 1 k + 2) • • • (n - 1 n)(nif) 



+ 



(4.34) 



One can straightforwardly obtain the higher terms, those terms higher in the number of 
propagators or UHV vertices. The total number of terms involved in the expansion (4.34) 
can be calculated as 



n—3 



E n-yPk — 2 



n—3 



(4.35) 



k=0 



where n-sCk denotes the number of /c-combinations out of (n — 3) elements which is also 
denoted as C(n - 3, k). As expected, this is equivalent to that of the expression (4.28) since 
we can easily count it as (1 + 1)™~ 3 . 

As mentioned in section 3, there are no apparent sums over permutations of number 
indices, or braid traces, in the definition of the massive holonomy operator. Such sums are 
already taken account of in the product of iterated integrals (3.11). The relevant sum is 
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given by the homogeneous sum (3.10) which, if explicit, produces (n — 1)! terms for the 
n-point UHV tree amplitude. In fact, this is what happens in the MHV tree amplitudes of 
gluons as well since the n-point MHV amplitude has (n — 1)! terms due to its braid trace or 
the sum over a e <S n -i- 

The UHV and the MHV amplitudes are different in structure, the physical configuration 
spaces are respectively given by = C n /S n ^ 2 and = C n /S n , respectively. Yet it is 
interesting to see that the above factor of 2 n ~ 3 arises from the absorption of the braid trace 
in a sort of compensating manner. It is also intriguing to compare the numbers involving 
terms for the UHV and the MHV amplitudes. These are given by by 2 n ~ 3 and (n — 1)!, 
respectively. The logarithm of these can be evaluated as (n — 3) In 2 < (n — l)[ln(n — 1) — 1] 
for a large n. 

From the expansion (4.34) one can easily visualize the expansion or the clusterization of 
the UHV tree amplitudes in terms of the UHV vertices connected by the massive propagators. 
This expansion is exactly what has been found in the derivation of Kiermaier's expression 
(4.28) in comparison with previously known results for the UHV tree amplitudes [7, 9, 11]. 
In the following, we interpret this UHV expansion in a functional language. 

Contraction of massive scalar operators 

As in the CSW rules for gluons, we can and should introduce a contraction operator 
involving the massive scalar operators. We can, for example, contract a pair of and a^°j 
(2 < J < n — 2) to replace it by a massive scalar propagator, with its momentum transfer 
given by pj in (4.29) or (4.30). Notice that once J is chosen, we consider the numbering 
indices in modulo J, i.e., 

J = -J = (mod J) (4.36) 
and generally J + i = — J + i = i (mod J) for < i < J. 

In analogy to the CSW rules (2.34), such a contraction operator can be defined as 



W {0) (x) 



cxp 



d/j,(Pj 



6af 8a ( _j 



,(o) 



,-iPj-(x-y) 



(4.37) 



V-*x 



where the limit y — > x is taken so that the time ordering x° > y° is preserved. The contraction 
operator can thus be expressed as 



W {0) (x) = exp 



d A P.i 



5 



-iPj-{x-y) 



(2n)*p*-mi \5af 5a^ J/ 
In deriving the above, we use the well-known identity 

d 4 q 



(4.38) 



dfi(q) e(x [) -y u )e 



_ vP^o-Mx-y) 



_ ^°\J<i( x -y) 



+ 9(y u -x")e 



-iq-(x-y) 



(2tt) 4 q 2 - m 2 + ie 

where q^ is an off-shell four-momentum with mass m and e is a positive infinitesimal. 
An S-matrix functional for the UHV tree amplitudes 



(4.39) 



21 



We now apply the Wick-like contraction operator (4.38) to the generating functional 
(4.23) for the UHV vertices: 



UHV 



= m°\x) J^ cx) \ a W c ,a<®] . (4.40) 



By construction, we then find that this functional serves as an S-matrix functional for the 
UHV tree amplitudes. Explicitly, the UHV amplitudes in the x-space representation are 
generated as 



a(±)c=a(°)=0 



^'"'-'"(x) - A§&\x), (4.41) 

/n— 1 

d/x(ft) II W^foMSS^W) , (4.42) 

A^ n \u,u) = -^- 2 (27r)V 4 )^ 1 + £p i+ p^ A^\u), (4.43) 

A^\u) = Tr(t c H c *---t c ^)C^\u), (4.44) 

° UHV {U) ~ (12)(23)-.-(n-ln)(nl) (nij " ^ UHV (°) ^ 



where we use the notation (4.32) in the last equation. Notice that here we are considering 
(nl) and ^u H v(o) 9n m ^ e ^-space representation. Thus all the massive propagators in 



(4.33) and (4.44) should be replaced by 

% d 4 P, 



Pj - m 2 (4tt) 4 P] - m 2 



(4.46) 



or simply by —idfi(Pj) in the expression (4.45). The extra factor —i here comes from the 
definition of the UHV vertex (4.25) in terms of the generating functional (4.23). Owing to 
the saturation of Grassmann variables, there arise no loop amplitudes made of the massive 
propagators for the UHV amplitudes 3 . Notice also that the UHV amplitudes preserve the 
single-trace structure due to the U(l) color degrees of freedom we assign to the massive 
scalars. 

Since the UHV S-matrix functional (4.40) leads to Kiermaier's expression for the UHV 
tree amplitudes, the above formulation shows nothing but a systematic functional derivation 



3 The massive loop structure can, however, enter in purely gluonic part of the UHV amplitudes. For 
example, we can easily consider massive one-loop subamplitudes for all-plus gluon legs; these are relevant to 
the so-called one-loop all-plus amplitudes in QCD. To incorporate these into the UHV amplitudes, we need 
to introduce massless propagators and an NUHV vertex. Thus such quantum effects do not arise as long as 
we use the UHV S-matrix functional (4.40). If we modify the S-matrix to include these ingredients as we 
shall do in the next section, however, the massive loop effects do arise in purely gluonic part of the UHV 
amplitudes. 
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of the massive CSW rules within the holonomy formalism, at least for the UHV amplitudes. 
We shall consider generalization to non-UHV cases in the next section. 

Lastly, we comment that in our formalism the mass effect does not break the supersym- 
metry. As analyzed in this section, the full M = 4 supersymmetry is crucial to derive the 
massive CSW rules. The form of the UHV vertex (4.14) suggests breaking of supersymmetry 
from Af = 4 down to Af = 2 if we treat the mass square in an isolated fashion. We have tried 
to implement such an interpretation, e.g., by considering a peculiar off-shell continuation of 
the Nair measure or by taking a different definition for supersymmetric massive operators, 
but any attempts did not work well. We thus come to realize that the holonomy formal- 
ism, or essentially Nair's prescription for superamplitudes, can naturally be continued to a 
massive system without breaking the extended Af = 4 supersymmetry. 



An S-matrix functional for non-UHV tree ampli- 
tudes 



As in the non-MHV amplitudes of gluons, the n-point non-UHV massive scalar amplitudes 
can be categorized as N fc UHV amplitudes in terms of the number of negative-helicity gluons, 
k — 1, 2, • • • , n — 3. The categorization is entirely gluonic since it dose not change the physical 
information of massive scalars but that of gluons. This means that such a categorization 
can be carried out by the original CSW rules or the purely gluonic MHV rules. Thus an 
S-matrix functional for massive scalar amplitudes in general can be constructed as 
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UHV 



[a^ c , a<°>] 



(h)c 



W^(x)W^\x)^ vcrtcx) 



UHV 



MHV 



(h)c 



(5.1) 



-^-(vertex) 
•> UHV 



[aW c ,a(°)] and J'mhv [a (h)c ] are given by (4.37), (2.34), 



where W (0 \x), W^ A \x) 
(4.23) and (2.33), respectively. 

It is tempting to construct the S-matrix functional without use of 7™ 

a O)d s i nce) as 

discussed in (3.6)-(3.8), the massive holonomy operator might include the purely gluonic 
holonomy structure. However, once we fix the physical configuration space for the 
massive scalar system, such an inclusion becomes physically difficult. For example, we need 
to separately define the braid traces of gluonic and massive part of the operators. There 
may be a way to circumvent these problems mathematically but it seems too artificial and 
so far we have not found any suitable methods that would lead to a definition better than 
(5.1). 

The NUHV tree amplitudes 

In what follows, we consider the next-to-UHV (NUHV) tree amplitudes, i.e., the simplest 
non-UHV amplitudes that contain a pair of massive scalars, one negative-helicity gluon, and 
an arbitrary number of positive-helicity gluons. Using the S-matrix functional (5.1), we 

can straightforwardly calculate the holomoprhic NUHV tree amplitudes ^NUHvro)( u )> the 
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counterpart of A^y^ (u) for the NUHV tree amplitudes, as 

n— In— 3 _ -i 

A(<t>i9aK) (,.\ _ 2"((»+'-+l)+-</ , n0i-(»-l)+'+)^ ,\ 1 2(H)- (»)+-o--(i+r)+)/ \ 

^NUHV^)^"^ ~~ ^UHVCO) I "7 2 ^MHVCO) I u ) 
i=2 r=l Hii+r 

n—1 n— 3 

= EE E^(* ffi -- , * <Ti+r * <+r+1 -- , * <_1 ) 

i=2 r=l <7G<S r +i 

P}((i+r+l)+-0„0i-(i-l) + i + )/ \ 1 (i)+-a_-(i+r)+)/ n 

°UHV(0) ^ a 2 MHV(O) \ u ,") 



Ul=U ii+r 



(5.2) 

ui=u a . a . A 



where, as usual, we consider the numbering indices in modulo n. The C's are therefore 
written as 



mi+r+i)+-4> n fa-(i-i)+i+) ( s m (nl)(nl) 

/UHV (°) W ' {i + r + li + r + 2)---{i-ll){li + r + l)' 



M(-l).(i) + ...a-...(i+r)+) ( x _ W , r n 

°mhv(o) i M > a ) — a„\f„„ \<„ „ \ i„ j\ y -^) 



where the off-shell momentum transfer and an associated spinor momentum are defined as 

iti+r = + Pf+l + ' ' ' + Vt+r = Pti+r + W7 lti+r 5 (5-5) 



pf&r = uf i+r u? i+r = ufuf. (5.6) 

Here w is a real number and r/^+ r is a reference null-vector. Permutation of the numbering 
indices for gluons is represented by 

a=( 1 1 + 1 l + T ). (5.7) 

\ &i ' ' ' O'i+r J 

Accordingly, the indices of the momentum transfer (5.5) are labeled by <7j's under the per- 
mutation a G iSy+i- In the expression (5.2), we denote ua +r = ui, or u CiCi+r = ui under the 
permutation, for simplicity. Also the SU(N) generators t c ^'s are abbreviated by t av s. 

Notice that the structure of the NUHV tree amplitudes is the same as the NMHV tree 
amplitudes of gluons except that in the former cases one of the MHV vertices is replaced by 
the UHV vertex. Consequently, there appear no sums over permutations over the numbering 
indices involving the UHV vertex. 

Examples of the NUHV tree amplitudes and comparison with BS expressions 

By construction, there are no 3-point NUHV tree amplitudes. Non-vanishing NUHV tree 
amplitudes start from n > 4. In what follows we consider first few examples of these. For 
n = 4, we can write down the NUHV tree amplitude as 

2(fa9t9a<M( \ £(04&J+)/ \ J_ 2(H)- 2+3_)/ x 



Q23 
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24 
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Tr(t 2 t 3 ) 
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2 °MHV(0) 1 
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+ 



(^3) 4 \ 



(l/)(/4)g 2 3 V(/2)(23)(3Z) (Z3)(32)(2Z); 



(5.8) 



Ui=U23 



where we use the fact that the color factor of the massive scalars is assigned to the U(l) 
direction of the U(N) = SU(N) x U(l) gauge group. Notice that the invariance under 
permutations of gluon legs is explicit in the above expression. In the literature this invariance 
is implicit and the amplitudes are usually given in a form of the first line in (5.8). Taking 
account of this fact, we find that the above expression agrees with the previously known 
result by Boels and Schwinn [15] with a certain choice of the reference spinor for fa. 

To be more precise, we can fix the reference spinor such that the on-shell partner pf of 
Pi is proportional to p% or || where p 2 denotes the gluon four- momentum with the 
numbering index 2. By doing so, we can easily see that the 4-point NUHV tree amplitudes in 
equation (2.10) of [15], the Boels-Schwinn (BS) expression, reduces to the above expression 
(5.8) since those terms proportional to {uiu 2 ) = (12) vanishes in the BS expression. Notice 
that the BS expression is stripped of color factors and permutation invariance under gluon 
transpositions. To compare the BS expression with our result (5.8), note also that the BS 
expression has a different helicity configuration from (5.8); the negative helicity is assigned 
to g 3 , not to g 2 , in there. 

The reference spinor of the above choice can be specified by the one satisfying rji = 
2c<yV ^2^ (Pi — CP2) where c is a constant. In the BS expression, the reference spinors for the 
massive scalars are set to identical, contrary to our choice in (4.20)-(4.22). Thus one can 
effectively reduce the BS expression to a form which is more compact than (5.8) by choosing 
a suitable reference spinor. In fact, such a choice was made in equation (2.11) of [15]. What 
we have shown here is that an alternative choice of the reference spinor leads to a different 
reduction of the BS expression where only the terms that involve the massless propagators 
survive. This reflects our basic relations in (4.18) and (4.19), i.e., the vanishing of the NUHV 
and UHV 2 vertices, in contrast to the BS relations in (4.15) and (4.16). 



d UHV 2 vertices, in contrast to the BS relations in (4.15) and (4. 
For n = 5 we can similarly compute the NUHV tree amplitudes 
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(5.9) 



M;=M 2 4 

(5.10) 



The rest of the 5-point NUHV amplitudes, A { ^^ 4 ^(u), can also be obtained from 
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symmetry argument. The expression (5.9) is in accord with the previously known result [16] 
with a certain choice of the reference spinor for (f>i in the same sense that we have argued in 
the case of 4-point amplitudes. 

To be more concrete, one can reduce the Boels-Schwinn expression for the 5-point NUHV 
tree amplitude, the one given in equation (3.10) of [16], to the form of (5.9) by choosing the 
on-shell partner of the massive scalar (pi to be proportional to p% in the BS expression. As 
in the case of the 4-point amplitude, such a choice removes all the contributions that do not 
contain massless propagators transferred by virtual gluons, and leads to the expression (5.9) 
once proper color structure and permutation invariance under gluon exchanges are imposed. 

As studied in [15, 16], the BS expressions for the 4-point and the 5-point NUHV tree 
amplitudes numerically agree with other set of NUHV tree amplitudes [7, 9] obtained by 
BCFW-type recursion methods. It is interesting to find that the above analyses show the 
connection between our formalism and the BS expressions at least for the cases of (5.8) and 
(5.9). For the other NUHV amplitudes, there exist no BS-type amplitudes in the literature. 
In this sense, generalization of the above analyses to n > 5 is not clear at this stage, however, 
we shall observe the appropriateness of our formalism for arbitrary n in the next subsection 
from a different perspective. 

For further understanding of our formalism we now present the 6-point NUHV tree 
amplitudes below: 
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(5.12) 



As far as the author notices, there exist no explicit expressions for the NUHV tree am- 
plitudes beyond n = 5 except the recent calculation of the so-called one-minus amplitudes 
in QCD at one-loop level, i.e., one-loop amplitudes of one negative-helicity gluon and an 
arbitrary number of positive-helicity gluons, with the internal lines being massive scalar 
propagators. The calculation is carried out by Elvang, Freedman and Kiermaier (EFK) in 
[22], applying the massive CSW rules for the one-loop amplitudes in QCD. The resultant 
one-minus one-loop amplitudes can easily be rendered into a form of the NUHV tree ampli- 
tudes and we can compare them with our general expression (5.2). For this purpose, we now 
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briefly review the EFK results for the one-loop one-minus amplitudes and their applications 
to the NUHV tree amplitudes for arbitrary n. 

Comparison with the EFK representation for the NUHV tree amplitudes 

The EFK calculation for the one-loop one-minus amplitudes is given in a form of in- 
tegrand which corresponds to the holomorphic amplitudes ^nuhv(o) 9n ^ n \u). One should 
also notice that the dependence on the reference spinors of (internal) massive scalars are kept 
explicit in the EFK calculation, i.e., the reference spinors are kept unspecified throughout 
the calculation. This is for the purpose of verifying the absence of spurious poles by use of 
the 77-independence of the amplitudes; see [22] for details. 




Figure 1: Diagrams contributing to the one- loop one-minus amplitudes of n gluons 



Following the notation in [22], the EFK result for the n-point one- loop one- minus ampli- 
tudes/integrands can be summarized as 
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where the terms in the right-hand side respectively correspond to the diagrams in Figure 1 
and are explicitly given by 
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In (5.15) the one-loop all-plus integrand Icsw + \a + 1, • • • , b, P) is defined as a trace over 
the UHV tree amplitudes: 
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(5.17) 



where the trace is taken over the two-component spinors and Tr' is defined as Tr'X = 
TrX — Trl to subtract the identity factor Trl = 2 from (5.17). In the expressions (5.14)- 
(5.17), the overall prefactor N p is introduced in order to make the pure gluonic amplitudes to 
those of QCD including (massless) fermion contributions in the loop. Neglecting the fermion 
contributions, we can set N p — 1. 

To obtain the NUHV tree amplitudes from these EFK results of the one-loop one-minus 
QCD amplitudes, we should bear in mind the following two things. 



1. One is the fact that in a one- loop (MHV) diagram we can make at least one leg on 
each side of the diagram be collinear to each other. This is due to the freedom we have 
in the choice of the reference spinors involving the loop propagators. This is also a 
useful lesson we have learned from the one-loop calculations of M = 4 super Yang-Mills 
theory in the holonomy formalism, see [28] for details. 

2. The other thing is, as mentioned earlier, that the reference spinors involved in the EFK 
results (5.13)-(5.17) are not specified. This has been convenient to study the freedom 
from spurious poles in the one-loop amplitudes [22]. However, for our purposes, i.e. 
to compare our formulation of the NUHV tree amplitudes (5.2) to the EFK results 
(5.13)-(5.17), we no longer need to keep the reference spinors arbitrary. We can fix 
them in a suitable way before reducing the EFK results to the NUHV tree amplitudes. 

From the first condition, we can easily find that the ring integrand i^ms' + '* vanishes upon 
the choice of u a || u a+ \, i.e., (aa + 1) = 0. The first condition also implies that diagram (Hi) 
in Figure 1 can be treated as a tadpole-like diagram. This means that the integrand IgprV 
can be considered as a UHV-MHV type amplitude upon the reduction to the NUHV tree 
amplitude by cutting the massive loop apart. The UHV part of the tree amplitude then 
has the minimum three legs composed of two massive scalars and one internal virtual gluon 
(but no pure gluons) . Similarly, we can reduce the subtree integrand IguUre^ ^° UHV-MHV 
type amplitudes upon cutting the massive loop apart. The UHV part of the reduced NUHV 
tree amplitudes have more than three legs (including an arbitrary number of gluons). These 
analyses show that the EFK representation for the NUHV tree amplitudes is given in terms 
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of UHV and MHV vertices connected by massive scalar propagators. This description agrees 
with our construction of the NUHV tree amplitudes (5.2). In fact, installing information of 
color factors and permutation invariance under gluon transpositions, we find that the above 
EFK representation of the NUHV tree amplitudes exactly agrees with our formulation (5.2). 

Generalization to non-UHV tree amplitudes 

Field theoretically it is straightforward to obtain non-UHV tree amplitudes out of the 
S- matrix functional for massive scalar amplitudes in (5.1). We simply apply a sequence of 
functional derivatives of interest to the S-matrix functional and evaluate the derivatives as 
in (4.41). Because of the contraction operators the computation is entirely based on the 
massless and massive CSW rules except that we make use of vanishing non-UHV vertices 
(4.18), (4.19) in the massive part, rather than the original proposal (4.15), (4.16). As in 
the cases of non-MHV gluon amplitudes, we can hence obtain non-UHV massive scalar 
amplitudes in terms of the UHV and the MHV vertices or what we previously call the "UHV 
rules." An explicit form of the non-UHV amplitudes are as tedious and complicated as that 
of the non-MHV amplitudes. From the S-matrix functional (5.1), however, we can easily 
obtain a succinct recursive expression for the non-UHV tree amplitudes: 



^NkUHVCO) ^ 
n— 1 n— 3 

<r . ( 5 - i§ ) 

i=2 r=l ln+r 



111 — Hi i+r 



where k — 1, 2, • • • , n — 3 and the meanings of qa +r and ui are the same as in (5.2). 
Alternatively, we can also express the N fc UHV massive scalar amplitudes as 
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n—l n— 3 I 

EV 2(( i + r + 1 )+-'l > ri<l>i-(i-^)+i+)(..\ 1 2((-0- W+-1fc--(*+0+)/,,\ 
^Nk-iUHVtO) ^ ' a 2 MHV(O) \ a ) 

i=2 r=l Hi i+r 



(5.19) 



where we make the negative-helicity indices implicit in ^^-luirvfof™^ 1 "'^ 1 ^ +l+ \ u )- The 
latter expression (5.19) is recursive in terms of the massive part of the massive amplitudes, 
while the former (5.18) is in terms of the gluonic part of them. The gluonic N fc MHV tree 
amplitudes in (5.18) can be constructed by the original massless CSW rules, while the massive 
N fc UHV tree amplitudes in (5.19) is constructed only from the N fe_1 UHV counterpart in an 
inductive way. In this sense the two expressions reveal the structure behind the CSW rules 
and the BCFW recursion methods and indicate their equivalence in the calculations of the 
massive scalar amplitudes at tree level. 

Lastly we should comment on loop effects to the massive scalar amplitudes. As noted 
earlier, massive loop effects arise in purely gluonic part of the massive scalar amplitudes. 
For example, we can incorporate one-loop all-plus gluon configurations into the gluonic part 
of the UHV amplitudes. There are also massless loop effects contributing to the general 
gluon amplitudes. Therefore, taking the expression (5.18) for instance, we can in principle 
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calculate the one-loop massive scalar amplitudes as 
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(5.20) 



From this expression we may obtain explicit forms of the one-loop massive scalar amplitudes 
but that task is beyond the scope of the present paper and we shall leave it to future works. 



6 Concluding remarks 

One of the main purposes of this paper is to investigate whether an off-shell continuation of 
Nair's superamplitude method can be applied to a massive model, particularly to amplitudes 
of an arbitrary number of gluons and a pair of massive scalars (or massive scalar amplitudes), 
in a framework of the holonomy formalism. In the present paper we have affirmed this 
proposition by making a specific choice of the reference spinors for the massive scalars, see 
(4.20). This allows us to obtain a functional derivation of the so-called ultra helicity violating 
(UHV) vertices (4.14) of the massive scalar amplitudes, eventually leading us to the S- matrix 
functional (5.1) for the massive scalar amplitudes in general. 

An essential ingredient of the S-matrix functional is given by the massive holonomy oper- 
ator @^(m) we have defined in section 3. This operator is relevant to the generation of the 
massive part of the massive scalar amplitudes. In particular, together with the contraction 
operator W^\x) in (4.37), this operator generates the UHV tree amplitudes (4.41)-(4.45) 
in a form that was previously reported by Kiermaier [21], using the so-called massive CSW 
rules of Boels and Schwinn [15, 16]. 

One of the interesting features in our formulation is that a careful analysis of the braid 
trace in the construction of O^(tt) leads to no apparent sums over permutations of gluons 
in the color structure of the UHV tree amplitudes, or the massive part of the massive scalar 
amplitudes in general. The number of terms involving the UHV tree amplitudes, however, 
does not drastically decrease from that of the MHV tree amplitudes of gluons. This is due 
to the fact that actions of the contraction operator W^\x) do not alter the gluon helicity 
configurations for the massive scalar amplitudes. We have briefly presented a quantitative 
analysis of this fact in (4.35) and below. 

The holonomy formalism is neither Lagrangian nor Hamiltonian formalism so that we 
do/can not introduce potentials for the incorporation of massive particles. Information of 
mass is embedded into physical operators such that helicity or polarization of the particles of 
interest is in accord with the definition of the helicity operator (2.31). In practice, this can be 
carried out by considering an off-shell continuation of Nair's superamplitude method. In the 
present paper we strictly follow this idea, with its concrete realization given in (4.6). Notice 
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that our approach is philosophically different from other approaches found in the literature. 
For example, the massive CSW rules [15, 16] are derived from the Lagrangian formalism 
[32]- [36]. The earlier approaches [4]- [11], on the other hand, focus more on the application 
of the CSW rules or the BCFW recursion relations to massive models and its usage rather 
than its derivation from first principles. The holonomy formalism is therefore qualitatively 
new in the studies of massive scalar amplitudes and possibly leads to a completely new mass 
generation mechanism. Obviously it is worth investigating how massive fermions and massive 
bosons will be incorporated into the same framework. We shall consider such extensions in 
a forthcoming paper. 

Lastly, we would like to emphasize that in our formalism mass effects do not invoke 
supersymmetry breaking. The full M = 4 supersymmetry has been crucial to derive the 
correct form of the UHV vertex (4.14) which forms a basic building block for the massive 
scalar amplitudes. (Such a construction is referred to as the "UHV rules" in the above text.) 
As explicitly shown in (4.23)-(4.26), the derivation is given in a conventional functional 
method. Thus the holonomy formalism can naturally be continued to a massive model 
without breaking the extended Af = 4 supersymmetry. It would be useful to take account of 
this fact in the construction of more realistic massive models in the framework of holonomy 
formalism or, more broadly, in the four-dimensional spinor-helicity formalism. 
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